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ABSTRACT 

We found a very likely counterpart to the recently discovered hard X-ray source IGR J20 18+4043 in the 
multi-wavelength observations of the source field. The source, originally discovered in the 20-40 keV band, 
is now confidently detected also in the 40-80 keV band, with a flux of (1.4 ±0.4) x 10~^^ ergs cm~^ s~^. 
A 5 ks Swift observation of the IGR J20 18+4043 field revealed a hard point-like source with the observed 
0.5-10 keV flux of 3.4!g;^x lO'^^ erg cm'^ s'^ (90% confidence level) at a = 20^8^38.^55, 6 = +40°4r00.''4 
(with a 4. ''2 uncertainty). The combined Swift-Integral spectrum can be described by an absorbed power- 
law model with photon index T = 1.3 ± 0.2 and A^h = ^-^-I'l ^ ^^^^ cm~^. In archival optical and infrared 
data we found a slightly extended and highly absorbed object at the Swift source position. There is also an 
extended VLA 1.4 GHz source peaked at a beam- width distance from the optical and X-ray positions. The 
observed morphology and multi wavelength spectra of IGR J20 18+4043 are consistent with those expected for 
an obscured accreting object, i.e. an AGN or a Galactic X-ray binary. The identification suggests possible 
connection of IGR J20 18+4043 to the bright 7-ray source GEV J2020+4023 detected by COS B and CGRO 
EGRET. 

Subject headings: ISM: individual (7-Cygni )— X-rays: individual (IGR J2018+4043; 3EG J2020+4017) 



1. INTRODUCTION 

The first hard X-ray imaging observation of the Cygnus 
region with the IBIS-ISGRI imager (Lebrun et al. 2003) 
aboard the International Gamma-Ray Astrophysics Labora- 
tory (INTEGRAL; Winkler et al. 2003) has revealed a new 
source, IGR J201 8+4043 (Bykov et al. 2004). The 320 ks 
fully coded field of view (FCFOV) observation yielded the 
source flux of (1.7 ±0.4) x 10~^^ ergs cm~^ s~^ in the 25- 
40 keV band. The Cygnus region is very rich in starform- 
ing sites. X-ray binaries, and supernova remnants (SNRs). 
IGR J201 8+4043 is located in the field of the SNR G78.2+2.1 
(commonly known as 7-Cygni SNR). A high-energy 7-ray 
source, 2CG 078+2, in the 7-Cygni field was discovered 
with the COS B satellite (e.g., Swanenburg et al. 1981). 
Observations with CGRO EGRET showed that this source 
(2EGJ2020+4026 = 3EG J2020+4017 = GEV J2020+4023), 
whose photon flux is ~ 1.2 x 10~^ cm~^ s~^ above 100 MeV 
(Sturner & Dermer 1995; Hartman et al. 1999), is the bright- 
est among apparently steady unidentified EGRET sources. 
Whipple 7-ray telescope observations (e.g., Buckley et al. 
1998) have established an upper limit of 2.2 x 10~^^ cm~^ s~^ 
for the flux above 300 GeV, indicating a break in the spec- 
trum above a few GeV. Bykov et al. (2004) suggested that the 
EGRET and ISGRI sources could be the 7-ray and hard X-ray 
counterparts of the same source in 7-Cygni , e.g., a Geminga- 
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like pulsar or extended emission from particles accelerated in 
the supernova ejecta/shell. In particular, the interaction of the 
wind of the early type 09V star HD 193322, located just 7' 
from IGR J20 18+4043, with 7-Cygni would be a plausible 
mechanism of particle acceleration. On the other hand, the 
ISGRI source position within the SNR field does not neces- 
sarily imply their physical connection. An example of such a 
confusion is the source IGR J17204-3554 projected onto the 
rich Galactic starforming region NGC 6334. The most likely 
identification of that source is a highly obscured extragalac- 
tic object NGC 6334B (see Bassani et al. 2005; Bykov et al. 
2006). 

Given the nominal 12' angular resolution (FWHM) and 
about V point source localization of ISGRI (Lebrun et al. 
2003), a reliable identification of the relatively weak (4-5 
a) ISGRI source requires high-resolution observations with 
other instruments. We report here the results of recent 1.5 
Ms FCFOV observations of the source with ISGRI (§2.1) 
and the detection of a counterpart of IGR J20 18+4043 in the 
0.5-10 keV band with Swift XRT (§2.2). The accurate po- 
sitioning of the Swift source allowed us to find its possible 
radio/IR/optical counterparts in archival data (§2.3 and §2.4). 
Based on these multiwavelength observations, we suggest that 
the source is a strongly obscured accreting source, likely an 
AGN, and discuss the possibility that the EGRET source is 
indeed its 7-ray counterpart (§3). 

2. OBSERVATIONS AND DATA ANALYSIS 

2.1. ISGRI observation 

The field of 7-Cygni has been frequently observed with 
INTEGRAL ISGRI. We used 1.5 Ms of ISGRI good time 
archival FCFOV data taken at revolutions 19-255 (2002 Dec 
9 - 2005 Nov 16). We applied the selection criteria taking 
only the science windows where the count rate was within 
±1 cr around the mean value determined by the whole 1 .5 Ms 
exposure. The standard red uction was pe rformed with the 
OS A 5.1 software (Goldwurm et alJl2003h . Pixel spreading 
was switched off for a more accurate flux estimation. Broad 
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Fig. 1 .— 40-80 keV ISGRI significance map (> 3.5 a) of the 7-Cygni field, with VLA 1.4 GHz contours overlaid. The Hght and dark 5^ x 5^ ISGRI pixels 
correspond to 3.6 and 3.8 cr significance, respectively. The VLA contours run from 9 to 73 mJy/beam on a linear scale. The Swift source position is marked by 
the white cross. The smoothed Swift XRT 2-8 keV image is shown in the inset. The bright 09V star HD 193322 was removed from the XRT image. 



energy bands were used to increase S/N. The source was de- 
tected with a flux of (1.2±0.2) xlO"^^ erg cm"^ s"^ in the 
20-40 keV band, and (1.4±0.4)x 10"^^ erg cm'^ s"^ in the 
40-80 keV band, obtained by cross caHbration with the Crab 
(see, e.g., Lubinski et al. 2004). The flux in 20- 40 keV band 
is below that quoted in Bykov et al. (2004) due to the new 
background model used in OSA 5 and more conservative se- 
lection criteria that we applied now to the much longer ex- 
posure dataset available. The 40-80 keV ISGRI significance 
map is presented in Figure [I] 

To search for ISGRI flux variations, we grouped all events 
in three time bins (the maximum number allowed by the count 
statistics), but found no significant flux variations. With the 
scarce statistics available, we can only exclude flux variations 
greater than a factor of 4 on a timescale of about a year. 

2.2. Swift observation 

The field of IGR J20 18+4043 was observed with Swift 
XRT (Burrows etal. 2000) on 2006 March 26, 29, and 30, 
with a total exposure of 4.8 ks ( Kennea et al. 2006). The XRT 
was operated in photon counting mode. The data were re- 
duced with the standard HEAsoft package (ver. 6.0.5). The 



event files were processed with the xselect tool (grades 0- 
12 were used). The XSPEC package (ver. 12.1.1) was used 
for the spectral analysis. A point-like source was detected in 
the XRT image at a position consistent with that of the IS- 
GRI source (see Fig. [D: a=20M 8^38.^55, (5=+40°4r00f'4 
(J2000), with a position uncertainty of 4'/2 (errors here and 
below in this subsection are at the 90% confidence level). 

We extracted 105 source counts from an aperture of 47 '.'2 
radius and estimated the observed 0.5-10 keV flux Fx = 
3.4;^Qg X 10~^^ erg cm~^ s~^. We fit the spectrum with the ab- 
sorbed power law model using the C-statistics and obtained 
the photon index F = l.l^os hydrogen column density 
A^H = 5.2!^;^ X 10^2 cm-2. The unabsorbed 0.5-10 keV flux 
is Fx = 5At_l^^ X 10-^2 erg cm'^ s'^. 

The combined Swift-ISGRI spectrum can be fitted with a 
single power-law model with F = 1.3!^;^, A^h = 6.1!^;^ x 10^^ 
cm~^. Based on the positional coincidence and spectral com- 
patibility, we conclude that we detected a Swift counterpart of 
IGR J20 18+4043; we wiU caU this source J2018 hereafter. 

Because of the presence of the bright (my = 5.82) 09V star 
HD193322 in the FOV, only W2 and M2 ultraviolet filters 
were used with the Swift/UYOT camera. The source was not 
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detected in these filters, the upper Hmits being not constrain- 
ing given the strong ISM extinction in this direction (see §3). 

2.3. VLA data analysis 

The radio image at 1.4 GHz was produced from VLA 
archival data obtained in D configuration on 1996 Sep 19. The 
data were processed with the MIRIAD software following 
standard procedures. The source 141 1+522 was used as a pri- 
mary flux density calibrator (5'i.4ghz = 22.8 Jy) and 1924+334 
as a phase calibrator. The synthesized beam is 54'' x 53'' 
(beamwidth at half power), P. A. = 33°, and the rms noise is 
about 5 mJy/beam. 

Figure 1 shows a contour image of the region around J201 8, 
at 1.4 GHz. The J2018 position is within an extended radio 
feature, of about 5' in size, ^ 70" south-east of the apparent 
position of radio intensity peak. The total flux density of the 
radio feature is about 480 mJy. The resolution in D config- 
uration is too low to resolve a possible point source within 
the structured extended feature. Overall, we cannot rule out 
the possibilty that the radio source is a counterpart of J2018, 
although the connection is not very certain yet. 

2.4. Optical and infrared data analysis 

Having the position of J2018 well constrained by the Swift 
observations, we looked for counterparts of the X-ray source 
in archival optical and IR data. We found a likely counterpart, 
with the flux of 88zb9 mJy at 8.28 /im, in the 6th MSX cat- 
alog (lEganetal .112001 at a = 20^8^39.4, 5 = +40°40'56'' 
(about 2" positional accuracy). We also found an extended 
MR source centered at a = 20^8^38.^73, S = +40°4r00f'l 
(2'.'6 from the Swiit position) in the 2MASS XSC catalog 
jSkrutskieetalJl2006l) . with magnitudes / = 13.0, if = 11.5, 
Ks = 10.7 (Fig. 2). These magnitudes are brighter than re- 
ported by Kennea et al. (2006) because those authors did not 
account for source extension. 

The optical images, taken from the POSS-II survey 
(lRddetal.lll991h . are also shown in Figure |2| A counter- 
part to J2018 can be seen in the POSS-Red (5878-7121 A) 
and POSS-IR (6939-9030 A) filters at the same position as 
the 2MASS source, and with a similar spatial extent. 

We estimated the fluxes of the extended optical source, us- 
ing the magnitudes of two nearby stars from the GSC 2.3 
catalog (iMc Lean et al.ir2000 ). marked with blue crosses in 
panels b and c of Figure 2. We estimated mRed =18.1 and 
^iR = 17.3 via aperture count-rate calculations. The absence 
of the source in the POSS-II Blue image yields an upper limit 
^Biue > 22.5. The aperture boundaries were set at the count 
rate corresponding to half a difference between the bright- 
est source pixel value and the local background level. Two 
reference stars in each band were used to estimate the zero- 
point and the count-rate vs. magnitude dependence. The 
same local background was subtracted from all the count 
rates, normalized by the aperture size. To convert the mag- 
nitudes into fluxes, we used the wab and POSS-II magni- 
tudes of the standard secondary spectrophotometric objects 
HZ 4 and G191B2B ( Oke 1990) and obtained Fgiue < 2.5/iJy, 
^Red = 0.27 mJy, Fir = 0.76 mJy for the optical source. We 
estimate the error of the cross-calibration, about 20%, to be 
mainly due to the extended nature of the source. 

3. DISCUSSION 

The Swift-ISGRI spectrum of J2018 shows that it is a non- 
thermal, highly obscured source. The hydrogen column den- 
sity estimated from the Swift data (§2.2) looks too high for an 



object in the 7-Cygni SNR (d ^ 1.5 kpc, expected A^h ^ 1- 
3 X 10^1 cm"2; Landecker et al. 1980). J2018 is likely a more 
distant object, not associated with 7-Cygni . 

A number of ISGRI sources have been identified as Galactic 
X-ray binaries, mostly HMXBs (e.g., FilHarte & Chaty 2004; 
Tomsick et al. 2006, and references therein). The J2018's X- 
ray spectrum (§2.2) and luminosity, ~ 3 x 10^^(J/10kpc)^ 
ergs s~^ in 1-80 keV range, as well as the magnitudes of 
the IR-optical counterpart (§2.4) do not contradict such an 
interpretation. The apparent ~10'' extension of the optical 
and 2MASS counterparts could be associated with a powerful 
wind of an HMXB's O or B companion. Alternatively, the ex- 
tended IR-optical source could be interpreted as a galaxy that 
hosts an AGN detected in X-rays. 

An extrapolation of the X-ray spectrum of J2018 into the 
EGRET range reaches the flux level of GEV J2020+4023, 
which suggests that these sources could be X-ray and 7- 
ray counterparts of the same object. An extensive dedicated 
search for a counterpart to 3EG J2020+4017 with Chandra 
and the Green Bank radio telescope in a close vicinity of its 
nominal position failed to show any clear candidate object 
(Becker et al. 2004; Weisskopf et al. 2006). However, J2018 
is located outside the fields of these searches. 

The position of J2018 is 30'.6 (35'.5) offset from the nomi- 
nal position of GEV J2020+4023 (3EG J2020+4017), exceed- 
ing the 8' (9 '.6) 95% localization accuracy in those catalogs 
(Lamb & Macomb 1997; Hartman et al. 1999). However, 
such offsets are not unusual among EGRET sources. The 
position accuracy of relatively faint EGRET sources in the 
Galactic plane was estimated by Gehrels et al. (2001) as ~1°. 
Moreover, even the very bright Vela and Geminga pulsars are 
located outside their EGRET 99% confidence contours (see 
Hartman et al. 1999). 3EG J2020+4017 is a bright source, 
but it is located in a very crowded region of the plane, with 
enhanced, highly inhomogeneous diffuse emission from the 
7-Cygni SNR and nearby molecular clouds, and at least two 
other EGRET sources (3EG J2021+3716 and J2016+3657) 
within the 3° vicinity, which may have affected the source 
localization. In such an environment, systematic uncertain- 
ties can increase the localization error up to ~0.°5. Thus, de- 
spite the positional discrepancy, we believe that one can con- 
sider the possibility that Swift-ISGRI and EGRET detected 
the same source. Future GLAST (e.g., Gehrels & Michelson 
1999) observation could resolve the issue. 

It is also worth mentioning the ^ 3a detection of a pos- 
sible Fe Kq, feature at 6.2 keV in the RXTE spectrum of 7- 
Cygni (Bykov et al. 2004). The feature, which could be at- 
tributed to a putative hidden accreting source (e.g., an AGN), 
does not appear in the ASCA spectra of the 7-Cygni regions 
observed by this satellite, which did not include the J2018 
position, while this position was well within the RXTE colli- 
mator field. 

To investigate the nature of J2018, we constructed its spec- 
tral energy distribution (SED; Fig. O using the ISGRI and 
Swift observations and the archival data described above. It 
should be noted that the data are noncontemporaneous, which 
is a point of some concern because the source could be vari- 
able. The Galactic absorption (Ay ^9.8) was accounted for 
in the unab sorbed NIR-o ptical spectrum using the reddening 
maps by Schlegel et al.' (1998)^. Being rather accurate for 
high latitude sources, the method often overestimates the red- 
dening by a factor of up to 1.5 for strongly reddenned objects 

^ See also |http://irsa.ipac.caltech.edu7applications/DUST/| 
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Fig. 2. — Archival POSS-II (top) and 2MASS (bottom) images of the J2018 field. The SwiftXKY position is indicated with the 4f^2 radius error circle. The 
green and blue crosses denote point-like sources from the GSC 2.3 catalog. 



(Arce & Goodman 1999). The optical emission is apparently 
extended, but with the data available it is not possible to disen- 
tangle the core component, corresponding to the point source 
detected by Swift, from the extended halo. Therefore, the 
optical points can be lower than shown in Figure 3, and the 
extinction corrected (with Ay = 9.8) optical fluxes in Figure 
13 should be regarded as upper limits, because of a possible 
overestimation (up to 3 mag) in the Ay used. 

The hypothesis that the 1 .4 GHz VLA excess is due to a 
point-like source at the Swift source position (e.g., an AGN) 
projected onto an extended radio source is not in conflict with 
the data available. Therefore, we conditionally include the 
VLA point, as well as the EGRET data, in Figure 3. Within a 
more conservative approach, they should be regarded as upper 
limits. 

The radio through hard X-ray SED in Figure|3lis consistent 
with that expected for a nearby Seyfert AGN or a radio-galaxy 
(e.g., Dermer & Gehrels 1995). For instance, at z = 0.02 the 
source luminosity below 100 keV would be about 10"^^ ergs 
s~^ 

If GEV J2020+4023 is indeed associated with J2018, then 
the radio through 7-ray source could be a blazar. Given the Ay 
uncertainties discussed above, the multiwavelength spectrum 
resembles the two-peak synchrotron self-Compton SED of S5 
0716+714, a blazar studied by Foschini et al. (2006). Si- 
multaneous multiwavelength observations are required, how- 
ever, to verify the blazar hypothesis and understand the pro- 



cesses responsible for such a spectrum. On the other hand, 
the leptonic jet model for a microquasar, such as discussed by 
Dermer & Bottcher (2006) for LS 5039, also provides a SED 
generally consistent with that given in Figure 3. 

NIR and optical spectroscopy could distinguish between the 
two interpretations by measuring the redshift of the galaxy 
that hosts the putative AGN. Also, microquasars are usually 
variable on minutes scale, while AGNs are variable on longer 
time scales making a variability study conclusive. 

The analysis given above allows us to conclude that the 
properties of the source discovered by Swift are consistent 
with those of IGR J20 18+4043 . However, that does not ex- 
clude a sizable extended hard X-ray flux from the possible in- 
teraction of the 7-Cygni SNR with stellar winds as discussed 
by Bykov et al. (2004). The extended VLA feature seen 
in Figure 1 can also be attributed to the SNR-wind interac- 
tion. High-resolution observations with Chandra and XMM- 
Newton, as well as simultaneous multiwavelength studies, are 
needed to separate the components and to make a firm conclu- 
sion on the nature of IGR J20 18+4043. 
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